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ABSTRACT 

We have analysed the Fermi LAT data on the SNR G73.9+0.9. We have confirmed a previous 
detection of high-energy y-rays from this source at a high significance of ^ 12cr. The observed 
spectrum shows a significant curvature, peaking in EFe at ~ 1 GeV. We have also calculated 
the flux upper limits in the mm-wavelength and X-ray ranges from Planck and XMM-Newton, 
respectively. We have inspected the intensity of the CO (1—>0) emission line and found a 
large peak at a velocity range corresponding to the previously estimated source distance of 
~4 kpc, which may indicate an association between a molecular cloud and the SNR. The 
y-ray emission appears due to interaction of accelerated particles within the SNR with the 
matter of the cloud. The most likely radiative process responsible for the y-ray emission is 
decay of neutral pions produced in ion-ion collisions. While a dominant leptonic origin of 
this emission can be ruled out, the relativistic electron population related to the observed radio 
flux will necessarily lead to a certain level of bremsstrahlung y-ray emission. Based on this 
broad-band modelling, we have developed a method to estimate the magnetic field, yielding 
B > 80/iG at our best estimate of the molecular cloud density (or less at a lower density). 
G73.9-I-0.9 appears similar, though somewhat weaker, to other SNRs interacting with a local 
dense medium detected by the LAT. 

Key words: acceleration of particles-gamma-rays: general-gamma-rays: ISM-ISM: super¬ 
nova remnants. 


1 INTRODUCTION 

G73.9-t-0.9 is a supernova remnant (SNR) located in the Cygnus 
arm of the Galaxy. It was first observed and id entified as an SNR 
in radio in the Effelsberg Galactic Plane sur vey i Reich_etaljl98^. 
High-resolution imaging at 1.42 GHz (see lKothes et al.ll2006l and 
references therein) shows a morphology consisting of a shell-like 
feature to the south with a sharp outer boundary and a centrally 
peaked diffuse emission to the east. 

The radio emission has the total flux of '-S Jy at 1.4 GHz 
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and the energy index of o- -0.23 (F^ oc v“; measured between 
0.4 and 10 GHz), corresponding to the photo n index of T = 1.23 
{F„ oc v‘ iLorimer. Lvne & Camilol (1 19981) searched for pulsed 
emission from the centre of the shell using the 76-m Lovell radio 
telescope, but without success. 

The region towards the source has been observed by th e 
Planck missiorQ jTauber et al.ll201(il : |Planck Collaborationll2015^ . 
but since G73.9-I-0.9 lies in a rather complex region of the Galactic 
plane, no results have been reported. Here, we obtain upper limits 
in the 30-143 GHz range, see Section l2n 

Optical line diffuse emission has been observed towards 


* wiki.cosmos.esa.int/planckpla2015/index.php/Main_Page 
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Figure 1. The colour image shows the map of the region in 0.3-3 GeV energy band with 0.1° x 0.1° pixels. All sources except 3FGL J2014.4+3606 
(G73.9+0.9) were included into the region model and thus are not present on the map; gives t he significance of a p oint-like source added to the model at 
a selected point. The green contours show the i mage of the SNR in 4.85 GHz continuum emission iCondon et alJl994i . The white contours are derived from 
the intensity data of the 115 GHz CO line (1^0; [D^ame et alj2Q0lh in the 1-7 km s * velocity range. 


the direction of this SNR llLozinskava. Sitnik & Pravdikovall 199.31 : 
iMavromatakisI l2003h . The optical spectrum implies an electron 
density of S 50 cm“^ at some locations within the SNR, and 
a moderate shock velocity of Vjh ~ 90 km s“'. These observations 
have allowed an estimate of the age of the SNR to be T ~ 11-12 
kyr. Its distanc e has been variously estimated. The most recent de¬ 
termination of IPavlovic et al.l l l2013ti puts it £> = 4.0 kpc, based on 
the radio surface-brightness to diamete r relation. The radial vel oc- 
ity of G73.9+0.9 has been measured bv iLozinskava et ^ ( Il993h as 
4 ± 3 km s“*. 

The SNR has not been detected in X-rays. We obtain an upper 
limit from XMM-Newton, see Section[2T| The region has also been 
observed with Imaging Cherenkov Telescopes at very high ener¬ 
gies, but no signal has been found. In parti cular, it was covered a s 
part of the HEGRA Galactic Plane survey jAharonian et alJl2002h . 
The 2-h exposure resulted in an upper limit on the number flux 
of N{>E) < 3.15 X 10“'^(£'/lTeV)“''^^cm“^s“‘ dAharonian et al.1 
l2002h . The VERITAS array also observed this region; the prelim¬ 
inary upper limit at the location o f the SNR is N{> 300 GeV) < 
2.3 X 10 cm-2 s ' ( lTheling|l2009li . 

At high-energy y-rays (0.1 GeV< E <100 GeV), 
iMalvshev, Zdziarski & Chemvakoval ( 1201 3l) have reported a detec¬ 
tion of y-ray emission from the direction of G73.9+0.9 using the 
Large Area Telescope (LAT) on board of Eermi. The y-ray emis¬ 
sion was discovered when investigating the region around Cyg 
X-1, which lies 3.2° away of G73. 9+0.9. The y-ray so urce was 
detected with the test statistic (Tj, iMattox et al.lll996l) of ==50, 


at RA 303.42, Dec. 36.21, which is only ~ 10' away from the 
catalogue position of the remnant, at RA 303.404 (20 13 37), 
Dec. 36.115 (36 06 54.0), corresponding to the Galactic coor¬ 
dinates of I = 73.77, b = 0.94. This detection has recently 
been confirmed in the Eermi Large Area Telescop e Third Source 
Catalogue (3FGL. Ipermi-LAT Collaborati^l201^ . as the source 
named J2014.4+3606 with the significanc^l of 4.5(T. Here, we re¬ 
analyse a significantly larger data set with the present instrumental 
calibration, and obtain significantly better statistical results. 

We present here a dedicated analysis of the currently available 
data set obtained with the LAT on the SNR and an investigation of 
the molecular content around it, derived from an analysis of the CO 
line (1—>0) intensity at 1 15 GHz (2.6 mm) in the region, observed 
with the CfA telescope dPame et alJI 19871) . We also study models 
possibly explaining the multi-wavelength spectrum of G73.9+0.9. 


2 MULTI-WAVELENGTH DATA ANALYSIS AND 
RESULTS 

2.1 The imaging and spectral analysis 

We have analysed the Fermi LAT data derived from a 15° radius 
region centred on the position of G73.9+0.9. Almost seven years of 
data obtained during MJD 54682-57209 (2008 Aug. 04-2015 Jul. 

^ heasarc.gsfc.nasa.gov/W3Browse/£ermi/£ermilpsc.html 
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06) were processed using the LAT standard tools (software version 
vlOrOpS), and analysed with the P8R2 response functions (SOURCE 
class photons). We have applied the standard cleaning to suppress 
the effect of the Earth albedo background, excluding time intervals 
when the Earth was in the field of view (EoV; when its centre was 
> 52° from the zenith), and those in which a part of the EoV was 
observed with the zenith angle > 90°. 

We have selected events with energies between 0.1 GeV to 
300 GeV, performing binned likelihood analysis. It is based on 
the fitting of a model of diffuse and point source emission within 
selected region to the data. The spatial model includes diffuse 
Galactic and extragalactic backgrounds (gll_iem_v06. fits and 
iso_P8R2_SOURCE_V6_v06.txt templates) and the sources from 
the 3EGL catalogue. We split the whole energy range into narrow 
energy bins, performing the fitting procedure in each bin separately. 
In each bin, the spectral shapes of all sources are assumed to be 
power laws with E = 2. The spectral shapes of diffuse Galactic 
and extragalactic backgrounds are given by the corresponding tem¬ 
plates. The normalisations of the fluxes of all sources and the dif¬ 
fuse background are treated as free parameters during the fitting. 
The analysis is performed with the python toolj^ The upper lim¬ 
its are calculated with the UpperLimits module provided with the 
FermilLAT software and correspond to the 95 per cent (=; 2.5cr) 
false-chance probability. 

We localise the position of G73.9-I-0.9 by building the Ts 
map of the selected region. The colour image in Pig. [T] shows 
this map, which illustrates the significance of adding a point 
source to the model of the region, w hich follows the dis¬ 
tribution and is approximately =; -fH jMattox et al.lll99al . The 
model in this case includes all the sources described above, ex¬ 
cept 3PGL J2014.4-t3606, which corresponds to G73.9+0.9 in the 
3FGL catalogue. We see a clear y-ray excess at the position of 
G73.9+0.9 in the 1-3 GeV energy band. The centroid of the im¬ 
age is the same as that in the 3FGL catalogue. The source shows a 
point-like morphology, with the 3cr upper limit on its radius (cor¬ 
responding to the decrease of Ts by 9 for the uniform-brightness 
disc model by 9 with respect to the point-source model) of 0.23° 
in the 0.3-3 GeV energy band. The maximum of the Ts for the 
y-ray source lies 0.1° from the centre of the radio emission from 
G73.9+0.9, which offset is not statistically significant. This is also 
shown in Fig.[T] where the black contours sh ow the 4.85 GHz con - 
tinuum emission from the GB6/PMN survey JCondon et alJll99^ . 

The energy spectrum for a point-like source centred at the po¬ 
sition of G73.9+0.9 was derived by means of the binned likelihood 
fitting, as described above. We show it in Fig.[2 The shown error 
bars correspond to Icr statistical errors. In the energy range of 0.1- 
300 GeV the spectrum can be fitted with a power law, which we 
show in Fig. 12 with the photon index of T = 2.73 ±0.13 and the 
normalization of No = (2.17 ± 1.10) x 10“** erg cm“^ s“*, defined 
as the value of EFe at 0.1 GeV. This fit yields the Ts ^ 64, which 
corresponds to cr =: 7.7. This confirms, at a hig her significance, 
the pre sence of this source in y-rays discovered bv lMalvshev et al.l 
( l2013h . 

However, the spectral points show a clear curvature in the 
spectrum, with the peak around 1 GeV. We confirm the high sta¬ 
tistical significance of the presence of a curvature by fitting a log¬ 
normal model (see AppendixlAt to the 0.1-300 GeV photon-energy 
band, which yields a much higher Ts - 150, which corresponds 
to the significance of ll.Str (with an addition of one free pa- 



Figure 2. The y-ray spectrum of a point-like source at the position of 
G73.9+0.9. The shown error bars present Icr statistical uncertainties and 
the upper limits correspond to 95 per cent confidence (about 2.5 (t). The 
dashed and dot-dashed curves shows the best-fit power-law and log-normal 
spectrum, respectively. 

rameter with respect to the power-law fit). The fit parameters are 
No = (7.7 ± 1.0) X 10“*^ cm“^ s“* MeV“* [defined as in equation 
(IA2 H.£h = 0.65 ± 0.04 GeV, which is the peak of the log-normal 
distribution in EFe, and p = 0.89 ± 0.07. We also obtain a simi¬ 
lar Ts for a broken power law fit, which has two more parameters 
than a single power law. We have also fitted the spectrum at £ > 1 
GeV by a power law, to facilitate comparison with Fer mi review pa- 

S iers, w hich give that index for other SNRs, see, e.g.. lBrandt et al.l 
201.5h . We have obtained T = 3.2 ± 0.2 and the normalization of 
No = (6.2 ± 0.8) X 10“*^ erg cm“^ s“*, given as EFe at 1 GeV. 

The integrated 1-100 GeV photon and energy fluxes derived 
directly from the obtained spectrum (i.e., treating the upper limits 
as zero points with the uncertainty treated as a systematic error) 
are == 2 . 0 +Q 24 x 10^^ cm“^ s“* and == 5.2 );q°j x 10“*^ erg cm“^ s“*, 
respectively. The estimated y-ray luminosities for the 0.1-100, 0.5- 
100 and 1-100 GeV ranges are, respectively, 28.3^5®, and 

9.9f\ X (D/4kpc)^ X 10^^ erg s“*. 

We have also searched for X-ray emission based on the XMM- 
Newton pointed observation of the field of G73.9+0.9 on 2006 Apr. 
19 (ObsID 0301880601, full window mode, PI: O. Kargaltsev). The 
observation yields an upper limit of £ 2.7 x 10“*^* erg cm“^ s“* 
in the energy range of 0.5^.5 keV within a circle of 90" radius 
around the observation pointing of RA 303.48, Dec. 36.25. The 
analysis has been done using the XMM-Newton online upper-limit 
softwarfl The upper limit is from the MOS2 detector, for which the 
exposure time is 2550 s; the MOSl and PN have much shorter ex¬ 
posures. From XMM-Newton and other X-ray observatories, there 
are no bright X-ray sources seen within the radius of 0.2° around 
the SNR centre. Fro m the image of the source in radio at 1.42 GHz 
jKothes et alj200^ . we have visually estimated the area of the peak 
of the SNR emission to be ~ (0.1°)^. This yields an upper limit on 
the SNR emission within 0.5^.5 keV of < 1.4 x 10“*^ erg cm“^ 
s“*. This corresponds to EFe < 6.2 x 10“*^* erg cm“^ s“* (assuming 
T = 2), shown in Fig. [2 This upper limit has a systematic uncer¬ 
tainty related to the uncertain source area and the spectral slope, 
which we roughly estimate to be of a factor of three. 

We have also searched for mm-wavelength emission using the 


^ fermi.gs£c.nasa.gov/ssc/data/analysis/scitools/python_tutorial itotanlLedas.ac.uk/flix/£lix.html 


© 2015 RAS, MNRAS 000,[T]{8] 























4 A. Zdziarski et al. 


Planck 2015 data dPlanck Collaboratioi]||2015bl O). The source does 
not appear in the Second Pla nck Catalogue of Compact Sources 
dPlanck Collahorationll2015^ . We have obtained 2cr upper limits 
in the 30-143 GHz frequency range using the 2015 full mission 
maps, and we show them in Fig.[3l 

The flux density of the source and the uncertainties used to 
calculate the upper limits have been measured as follows. For the 
flux density, aperture photometry centred at the position of the 
source has been used using an aperture with a radius equal to one 
FWHM. The flux densities have been corrected to take into account 
that a fraction of the beam solid angle falls outside the aperture. 
For the rms, we have used an annulus with the inner radius of 1 
FWHM and the outer radius of 2 FWHM. For all the frequencies, 
the FWHMs are the effect i ve beam s as listed in column 3 of Table 
2 of IPlanck CollaborationI d2015dh . 


2.2 A molecular cloud around G73.9+0.9 


Several GeV-bright SNRs have been shown to interact with co¬ 
located molecular clouds. The presence of such a cloud appears 
to be requi red for emission o f y-rays detectable by the LAT from 
old SNRs l lHewittetalJl2013h . A safe identification of interaction 
with such high d ensity material wo uld be a detection of masers 
(see, e.g., W51C. lOreen et aDl2006l) . A search in the direction of 
G73.9-I-0.9 for both OH and H 2 O masers, at 1.72 GHz and 22 
GHz, respectively, gav e no results faewitt & Yusef-Zadehll200^ : 
IWoodall & Gravl2007h . Data from the I-GALFA Hi 21 -cm line sur¬ 
vey did not show any si gnificant excess e mission that can be asso¬ 
ciated with the remnant jPark et al.ll2013h . 


Therefore, in order to evaluate the molecular content at the lo¬ 
cation of the remnant, we use the intensity, Hco, of the CO line 
(1—>0) at 115 GHz (2.6 mm), given in units of K km s ‘. CO is 
the most commonly used tracer for molecular material, given the 
lack of a possibility to directly observe H 2 . In order to obtain Ahj, 
the column density of H 2 , a co nversion factor Xqo is then used 
( iBolatto. Wolfire & Lerovll2013h. An extensive archive is provide d 


s p rovide d 

ai l l200lh . 


by iDame et al.l (19870 and iDame, Hartmann & ThaddeusI 
The data are releaseqj in cubes of radial velocity (in the local stan¬ 
dard of rest) and Galactic latitude and longitude. We can estimate 
the mass of the molecular material from Wco as 


M = ^D^AGp.Aco^Wco cx Am, 

px 


( 1 ) 


where X is the H mass fraction, nia is the mass of the H nucleon, 
and Xco is a conversion factor (see below), and Aflpx corresponds 
to the solid angle subtended for each pixel in the map (square bin¬ 
ning of 0.125° per side). The term 2px Wco takes into account the 
binning in velocity of the data cube and is obtained by summing 
the map content for the pixels in the desired sky region and desired 
velocity range and scaled by the bin size in velocity. 

We have calculated the absolute v alue of the molecul ar mass 
using the conversion factor estimated bv IPame et al.l (200 ih . Xco = 
1.8 X 10^° cm“^(Kkm/s)“'. Studies of the estimation of the gas 
content using differe nt tracers show compatible average results 
(Pedaletti et al.ll2015l) . 

We have found that the velocity distribution of the CO con¬ 
tent shows a clear excess in the range of velocity of (-2.0, 3.2) km 
s“*. This o verlaps with the v e locity range measurement of v == 1-7 
km s“* of iLozinskava et ^ (l993h . Thus, we adopt the range of 


^ www.cfa.harvard.edu/rtdc/CO/ 


V 1.0-3.2 km s“' as the SNR systemic velocity. Ba sed on it 
and us ing the Galacto-centric rotation-distance relation of iciemenj 
lll98.5h . we estimate the distance of G73.9+0.9 of either 0.2-0.5 kpc 
or4.3^.5 kpc. These estimates bear some sys tematic uncertai nties 
due to the uncertainty of the rotation curve of lClemeii sl( ll985t) and 
to the possible proper motion of the source. The lower range is 
highly unlikely, while the upper one approx imately agrees with the 
estimate of 4.0 kpc of IPavlovic et all ( l2013h . 

We have then estimated the mass content in a region of 0.15° 
sky integration radius around the position of the y-ray source as 
M 4300(D/4kpc)^MQ. The conversion of the mass to density is 
uncertain, given the lack of knowledge about the shape of the cloud. 
We provide a rough estimate of the electron density of the medium 
assuming a uniform spherical cloud with the radius equal to the sky 
integration radius, which is, obviously, a rather strong assumption. 
For X - 0.7, this yields 30(D/4kpc)~‘ cm~^. T his agrees with 
the upper limit of 50 cm“^ of iMavromatakisI (20031) . However, the 
actual density can be lower if the radius of the cloud is higher than 
that corresponding to 0.15°. Furthermore, we cannot exclude that 
the directional coincidence of the molecular cloud and the SNR is 
accidental. 


2.3 The hadronic model for the y-ray spectrum 


Given the relatively high density of the medium where the SNR 
is located, two standard models can be considered to explain the 
observed y-ray emission. In both, the emission is due to radia¬ 
tion of the particles accelerated in the supernova shock interact¬ 
ing with the dense medium around it. In the first model, the y- 
rays are produced through the decay of neutral pions produced 
in the collisions between the accelerated ions at high energies 
and the ambient ions. In the second model, the y-rays have a 
leptonic origin from non-thermal bremsstrahlung radiation of the 
accelerated non-thermal electrons on the ambient medium (see 
Section Imj. A possible contribution from inverse Compton ra¬ 
diation can be neglected at £ 100 GeV in comparison with the 
bremsstra hlung given the high ambient density we ha ve measured 
(see, e.g.. lYang. de Ona Wilhelmi & Aharonianll2014h . 

We es timate here the hadron ic y-ray emission using the for¬ 
malism in iKafexhiu et al.l (20141) . The observed y-ray spectrum 
is rather narrow, and we have found that a correspondingly nar¬ 
row proton distribution is required to reproduce it. Given that 
we have only four data points, we have not attempted a formal 
fit, but only tested a range of distribution. The observed spec¬ 
trum can be reproduced by a proton distribution with a power- 
law index of p = 2.0 and an exponential cutoff at 20 GeV, 
i.e., N(E) oc exp(-£/20GeV). Decay of 7r° produced in 
proton-proton collisions yields the spectrum shown by the solid 
curve in Fig. The energy content in the proton distribution 
above the threshold for pion production, i.e., E > 0.3 TeV, is 
4.7(D/4kpc)^(ne/30cm)“' x 10"^* erg. We note this is a mod¬ 
est energy content when comparing with the total kinetic energy of 
an SNR, even one with the age ^10 kyr. 

The steep high-energy part of the y-ray spectrum is sim¬ 
ilar to those found in a number of GeV-bright SNRs and can 
be explained naturally by either a deficit of high-energy protons 
due to the cooling or escape of particles during inefficient cos- 
mic ray acceleration in dens e surrounding environments (see, e.g. , 
Aharonian & AtovanI Il996l: iMalkov, Diarnond & Sagdee"^ 1201 iL 


2 OI 2 I : Gabici. Aharonian & Casanovtj 200^ . Alternatively, it can 
be a result of the spectral shape emerging from re-acceleration 
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Figure 3. The obsei'ved broad-band spectrum of G73.9+0.9. The black sy mbols with open cir cles give our y-ray measurements and upper limits, and the 
magenta symbols give our upper limits from Planck. The radio data from iKothes et alj j2006h and the upper limit from XMM-Newton are shown by blue 
crosses and error bars, and by the blue upper limit with an open circle, respectively. The spectrum from decay of neutral pions in the hadronic model is shown 
by the black solid curve, see Section lT^ The cyan dot-dashed curve shows the bremsstrahlung emission of e* from decay of charged pions in the same model 
assuming the SNR age of 12 kyr and at = 30 cm“^ and B = 80//G. The corresponding synchrotron emission is below the shown range of fluxes. The radio 
spectrum is modelled by synchrotron emission (shown by the red dotted curve) from e-folded power law electrons with the minimum value of By^ allowed by 
the radio data not showing any high-energy cutoff. The normalization of the spectrum yields liAnD^). The red dashed curve shows the maximum 

bremsstrahlung emission from the same electrons allowed by the data, which implies that B ^ 80//G. See Section [T4l for details. 


of pre-existing cosmic rays in co mpressed clouds jUchivama et al.l 
l20ld:lBlandford & Cowidll982h . 


We then consider the possible role of e* produced from decay 
of charged pions. We have calculated the rate of their production, 
Qe±{y), which is relatively similar both in shape and normalization 
to the photon production rate from Td decay. The rate Qc+Ct) con¬ 
tains most of the energy around the Lorentz factors of 7 ~ (100- 
2000). The e* lose then their energy in interactions with the back¬ 
ground gas via Coulomb and bremsstrahlung, and with the mag¬ 
netic field via the synchrotron process. We have calculated the en¬ 
ergy loss rates, 7 , for those processes assuming = 30 cm“^ and 
B = 80;rG (see Section l2!4t . The time scale for the total energy 
loss, yly, is found to be > 10^ yr for 7 =: (10^-5) x 10'*, which 
is the range relevant for the e± emission. This is much longer than 
the age of the SNR of T =: 11-12 kyr. Thus, the energy losses 
can be neglected and an upper limit on the current electron dis¬ 
tribution within the source can be estimated as N{y) — TQs±(y) 
at r = 12 kyr. This is an upper limit since we neglected possi¬ 
ble leakage of e* out of the source. We have then calculated the 
synchrotron and bremsstrahlung spectra for this ^( 7 ) (see Section 
l2.4l for the method). Both processes are found negligible. The latter 
component is shown by the dot-dashed curve in Fig. [3 whereas the 
former is below the scale of the plot. 


Our present modelling predicts r ather few photons at £ > 50 
GeV. This can be tested by the CTA dAcero et alj|2013l) . the next 
generation of Cherenkov telescopes. If there is still emission at this 
photon energy range, the high angular resolution of the CTA would 
permit a detailed investigation of the acceleration region. 


2.4 Constraints on the electron distribution and magnetic 
field from the broad-band spectrum 


As we have found out above, the proton distribution required to 
explain the 7 -ray spectrum is quite narrow. On the other hand, the 
observed lack of any visible curvature in the radio spectrum, see 
the black data points in Fig.O requires the existence of an electron 
distribution with a power law form for at least two decades, though 
the 70-GHz Planck upper limit requires some steepening by this 
frequency. This implies the protons and electrons are accelerated 
differently. 

The local electron power-law index of p ^ 2a + \ ^ 1.46 
is required. Then, the observational data impose constraints on the 
overall electron distribution and the magnetic field via the normal¬ 
ization of the radio spectrum (emitted via synchrotron), the lack of 
its curvature up to 10 GHz, and the 7 -ray detection and upper lim¬ 
its (considering the contribution from bremsstrahlung by the same 
electrons). 

We first assume the distribution of relativistic electrons in the 
source to be a generalized e-folded power law, 

N(y) = exp [-( 7 / 70 )^]. ( 2 ) 


We have found that either q = I (exponential cutoff) or q = 2 
(super-exponential cutoff) lead to almost the same constraints on 
the magnetic field. Thus, we present the results only for ^ = 1. We 
also consider electrons with a broken power-law shape, having the 
index of p + Ap above 70 . We assume a smooth transition between 
the two power-law parts. 


V(7) = 


yP -|- yf^yPP^P 


(3) 


Our procedure is as follows. The synchrotron power-law spec¬ 
trum constrains its index and a combination of the electron normal¬ 
ization and the magnetic field strength, B. The presence or absence 
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of a curvature constrains a combination of B and 70 . The resulting 
three parameters are 


P. 


^^g(l+p)/2 ^ ^ 

4;r£»2 " ‘ 


Byl ^ A,_. 


(4) 


The overall synchrotron spectrum depends only on these three pa¬ 
rameters, with its shape given by p and A 2 , and the normalization, 
by Al. The lack of any visible spectral curvature up to 10 GHz 
yields a lower limit of A 2 . A further constraint can be provided by 
the bremsstrahlung spectrum emitted by the same electrons, given 
the density of the background matter and the observational con¬ 
straints in y-rays. This allows us to constrain all the four individual 
parameters, p, B, and 70 . 

Although the local power-law index of the electrons responsi¬ 
ble for the observed radio emission has p ^ 1.46 and the normaliza¬ 
tion of Al 1.3 (cgs), a harder distribution is required in the pres¬ 
ence of an exponential cutoff, equation ([ 2 l( with ^ = 1 , which cutoff 
softens the local spectrum. We constrain the considered range of the 
index to p > 1 . 2 , given the comp lete lack of harder distributions in 
radio SNRs dKothes et al.ll200^ . At p = 1.2 and Ai 0.7 (cgs), 
A 2 - 5 X lO"* G corresponds to the (approximate) maximum spec¬ 
tral curvature required by the data at < 10 GHz and the minimum 
one required by the 70-GHz 2cr upper limit. 

We show the synchrotron spectrum corresponding to those 
constraints (with the minimum of By^) by the red dotted curve Fig. 
[3 Since we calculate the emission from the entire SNR, we use the 
for mula for the synchrotron emis s ivity averaged over the pitch an¬ 

gle dCrusius & Schlickeise3ll986l : lGhisellini. Guilbert & SvenssonI 
ligsaT which we then integrate over the electron distribution. For 

a value of A 2 higher than the above lower limit, the synchrotron 
spectrum in the observed radio range will be flatter, and the peak in 
vFy will be higher and at a higher energy than that plotted. 

We then calculate the bremsstrahlung spectrum from the same 
electrons, using our estimated density of the molecular cloud, 
Ue - 30cm“^. We use the bremsstrahlung formulae for £ 2; 0.1 
GeV, at which energie s the contributions from ions and ele ctrons 
have the same form dStrong, Moskalenko & Reimeill20o3) . The 
bremsstrahlung spectrum is oc A'e«e(3 + X)(l -t- Xy^{AnD-y^. Its 
photon index at £ <K yom^A is == p, and the peak of its EFe 
spectrum is around £ ~ m^cryg at roughly EF^ SS n^K^y^^'’. 
We then compare the calculated bremsstrahlung spectrum to the 
constraints from the Fermi observations. We find the 7 -ray ob¬ 
servations constrain B > 80 (D/ 4 kpc)“^/T+i)yjQ 

KJiAnD^) - A,/B<‘+'’W 2 ^ (A 2 /S)‘/" - 2.5 x lOl 

The model with Smin and 70 ,min is plotted in Fig.[3as the red dashed 
curve. We see the predicted bremsstrahlung component is at the 
observational upper limit. We have also found that if we decrease 
(increase) 70 via a decrease (increase) of A 2 , we need to harden 
(soften) the electron distribution in order to have the synchrotron 
spectrum passing through the data points. While the decrease of A 2 
decreases the resulting bremsstrahlung spectrum, the decrease of p 
increases it, and the value of B^^n remains almost unchanged. Thus, 
our constraint is robust, independent of details of the estimate of the 
possible cutoff in the synchrotron spectrum. 

We then repeated the calculations for electrons with a broken 
power-law distribution, equation 0 - However, we have found that 
unless Ap ^ 1.5, the resulting constraints impose an even higher 
Smin- At Ap =1.5 (and = 30cm“^), - 80 /jG, i.e., the same 

as for the exponential cutoffs. Thus, our constraint on B is again 
robust, approximately independent of the assumed break or cutoff 
in the electron distribution. 

However, the actual matter density within the SNR can be 


lower than 30cm“^, e.g., due to the cloud length along the line 
of sight being more than we assumed. Then, the bremsstrahlung 
emission will be weaker, and the minimum allowed magnetic field 
strength will be lower. A relatively strict lower limit on the density 
is provided by that of the ISM, which we take as Ue = 1 cm“^. In 
this case and for equation 0 with q = 1, we find Bmin - 9 pG and 
70 ,mm - 7.5 X 10“*. 

We point out that it would be of significant interest to detect 
G73.9-t0.9 at a frequency higher than the current detection at 10 
GHz, to improve the constraint on the electron distribution obtained 
in this paper. The Planck upper limit at 70 GHz already implies 
some steepening. 

We note that our results also imply that the observed 7 -ray 
spectrum cannot be due to the bremsstrahlung emission from the 
population of electrons emitting the radio spectrum. In order to 
increase the normalization of the relativistic electrons to the EFe 
level corresponding to the observed 7 -ray peak, we have to strongly 
decrease the magnetic field strength, to <K Smin- This then would 
imply m^c^yo » 1 GeV, i.e., the bremsstrahlung peak much above 
that observed. Alternatively, if we use a value of 70 correspond¬ 
ing to the observed 7 -ray peak, the synchrotron spectrum would 
break at an energy much below that of the highest point of the ra¬ 
dio power-law. 

Substantial work al ong the above lines has been 


done before. 


see, e.g., 


Abdo et al 


J2OO9L l2010tJl3) 


_^. . and 

ISrandt & Fermi-LAT CollaborationI ( Eoiik ^However, those 
authors have only obtained some specific sets of the parameters 
resulting in good fits to their data and did not apply the formalism 
outlined here. Our new method can be applied to other SNRs 
for which there is a good measurement of the radio spectrum in 
a wide frequency range, an estimate of the background matter 
density, and either a detection or upper limits in the 7 -ray range. 
The presented formalism uses the bremsstrahlung emission of the 
synchrotron-emitting electrons. However, the method can be easily 
generalized to include inverse-Compton emission of the same 
electrons, provided the seed photons for scattering are specified. 


2.5 Comparison to other sources 

We c ompare G73.9-H0.9 to other SNRs detected by the Ferm i 
LAT llThornn son. Baldini & Uchiya^ I2OI2I : iHewitt et alj l2013l : 
ISrandt et^IT 2015h . iHewitt et al.l ^201; ) find the 7 -ray SNRs fall 
into two categories: the young ones and ones interacting with a 
local dense medium (with n ~ 10^ cm“^). G73.94-0.9 appears 
to belo ng to the latter cate gory. On the radio vs. 7 -ray flux dia¬ 
gram of IHewitt et al.l ( l2013h . G73.94-0.9 falls below all other iden¬ 
tified interacting SNRs. However, this may be due to its rela¬ 
tively lar ge distance; the rang e of the 1-100 GeV luminosities also 
shown in lHewitt et al.l ^20131) is = (0.7-40) x 10^^ erg s“*, and for 
D = 4kpc, the £(1-100 GeV) of G73.9-^0.9 is in the middle of 
their range. Still, most of the sources shown in lHewitt et al.l l l2013l) 
have L above the middle value. The relatively low luminosity of 
G73.9-f0.9 may then be due to the relatively low density of its 
molecular medium. _ 

From iKothes et alj j2006l) . we find the SNR radio indices to be 
between a = -0.2 and -0.6, with the index of G73.9-I-0.9 being on 
the hard side. On the other hand, the 7 -ray index at £ > 1 GeV of 
= 3.2 ± 0.2 is on the soft side of the observed range of F 1.5-3.5, 
and overlapping with the range of F = 2.1-3.1 foun d for SNRs in¬ 
teracting with molecular clouds l lBrandt et alj201^. The pr ojected 
size of the SNR, with the 27' total diameter ( Greenll2oT4l) . corre¬ 


sponds to = 31(I)/4kpc) * pc. The LAT-detected interacting SNRs 
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have diameters of ==10-50 pc jThompson et aljUoi2l : lHewitt et alj 
l2013h . This comparison also favours the distance of D ~ 4 kpc for 
G73.9+0.9. 

Finally, we note that a number of pulsars (some of them 
mill isecond ones) hav e spectra similar to that of G73.9+0.9, 
see lAbdo et al.l ll2013h and the associated online spectr^E As 
some examples, we list PSR J0023+0923, PSR JOlOl-6422, PSR 
J0357+3205, PSR J0622+3749, PSR 10908^913, PSR J1057- 
5226, PSR J1744-1134, PSR J1833-1034, PSR J2030+3641. 
However, there is no evidence for the presence of a pulsar in 
G73.9+0.9, and the similarity may be incidental. 


3 CONCLUSIONS 

We have analysed almost seven years of Fermi LAT data on 
G73.9+0.9. We have confirmed the detectio n of high-energy y-rays 
from this source by iMalvshev et al.l (l2013h . but at a much higher 
statistical significance of = 12cr. The detected spectrum covers the 
0.1-5 GeV energy range and shows a significant spectral curvature, 
peaking in EFe at ~1 GeV. 

We have found that G73.9+0.9 is most likely located within 
a molecula r cloud, which distan ce we have estimated as ~4 kpc, 
confirming IPavlovic et al.l (l2013h . We have estimated its density as 
a function of the source distance. The y-ray spectrum can be fit¬ 
ted by emission from decay of pions produced in collisions of ions 
accelerated to relativistic energies with a power-law distribution in¬ 
teracting with the matter of the molecular cloud. 

While the leptonic origin of the y-ray emission can be ruled 
out, some bremsstrahlung y-rays have to be emitted by the relativis¬ 
tic electrons emitting the observed radio flux (via the synchrotron 
process). Using this, we have developed a method of estimating 
the magnetic field in the synchrotron-emitting region. It implies a 
lower limit of ~ 80yuG if the molecular cloud electron den¬ 
sity is =:30cm“^ (estimated by us assuming the molecular cloud 
is spherical with the radius corresponding to 0.15°), or less for 
a lower density. The magnetic-field estimate is robust for a given 
background density, almost independent of the shape of the radio 
spectrum above the measured range. 

We have compared G73.9+0.9 with other SNRs interacting 
with surrounding media detected by the LAT, and found the source 
to be relative similar though with weaker radio and y-ray emission. 


NOTE ADDED IN PROOE 


all three frequencies), and there is neither an optical nor an X-ray 
counterpart seen. Thus, the source does not appear to be a blazar. 
If the source is a typical FR II radio galaxy, its relatively small an¬ 
gular size would suggest it is much more dista nt than most of radio 
galaxies detected by Fermi jAbdo et al. l2010bh . We conclude that it 
is unlikely (though we cannot rule it out) that the source contributes 
noticeably to the observed GeV emission. 
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10-41og5o(£/£b)^ 


(A2) 


A similar form of the log-normal dis tribution was used to fit blazar 
spectra b y, e.g., Lan^^et^ (Il98fih . The values of Ao and £b are 
given by jMassaro et al.ll2003l . 


No = Eb = 


(A3) 


Their uncertainties then need to be calculated with the propagation 
of the errors on a and fj. On the other hand, setting a = 2 and 
leaving the other three parameters free converts the spectrum of 
equation dAlt to the log-normal spectrum in EF^, equation ( IA2b . 
In this case, Eo does give the peak of EFe ■ We use this form in our 
spectral fitting. 
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APPENDIX A: THE LOG-NORMAL/LOG-PARABOLIC 
SPECTRUM 

The Fermi EAT analysis softwar^^ provides a log-parabolic photon- 
spectrum model for spectral fitting following the form given by 
iMassaro et al.l ( |2004|) . 


dE “U' 


-a-/31og,(,(£/£') 


(Al) 


with the four above parameters. However, a parabola has only three 
independent parameters. In our case, an observed spectmm in the 
logarithmic space, i.e., logiQcLA/dE, can determine the three coef¬ 
ficients of log'iQ E, i = 0, 1 2. The square term has the coefficient 
of p, and thus it can be uniquely fitted given an observed spectrum. 
However, the constant and linear terms are given by combinations 
of all four variables in equation dAlb . Thus, N'g, a, and E^ cannot 
be separately determ ined, due to their internal correlation. Indeed, 
iMassaro et al.l ( |2004|) set EJ, fixed at 1 keV. 

However, if we keep fixed to some value (as recommended 
now on the Fermi web page), then it will not have the meaning of 


^ fermi.gs£c.nasa.gov/ssc/data/analysis/scitools/source_models.html 
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